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1 Introduction 

In the Standard Model (SM), a weak decay such as B° — > vV(+'y) can only occur 
through second-order diagrams like those shown in Fig. [TJ All these processes are 
highly suppressed within the SM. Like all purely leptonic B decays, they contain 
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Figure 1: The lowest-order SM Feynman graphs for B° — > invisible (+ 7) decays: a) 
box diagram, b) and c) weak annihilation diagrams. 

a b — > d transition plus an internal quark annihilation that further suppresses the 
amplitude with respect to rare semileptonic decays. In addition, helicity suppression 
factors proportional to ml make the uV channel completely undetectable in the SM 
scenario. For the 1/F7 channel the latter factor is not present, resulting in SM branch- 
ing fraction expectations at the 10 -10 level pp. Several new physics models predict 
enhancements of these branching ratios up to values close to the experimental detec- 
tion: a phenomenological model allows for associated neutralino production from B° 
decays with a branching fraction in the 10~ 7 to 10~ 6 range |2|. Also, models with 
large extra dimensions can have the effect of producing significant, although small, 
rates for invisible B° decays [SI IH E] • 

The data used in this analysis were collected with the BABAR detector at the 
PEP-II e + e~ collider at SLAC. The data sample corresponds to a luminosity of 
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424 fb _1 accumulated at the T(4S) resonance and contains (471 ± 3) x 10 6 BB 
pair events. For background studies we also used 45 fb _1 collected at a center-of- 
mass (CM) energy about 40 MeV below the BB threshold (off peak). A detailed 
description of the BABAR detector is presented in Ref. [5J. 
More details on these analysis may be found in Ref. [7]. 

2 Reconstruction and selection 

The detection of invisible B decays uses the fact that B mesons are created in pairs, 
due to flavor conservation in e + e~ interactions. We reconstruct events in which a 5° 
decays to D^ + i~T> (referred to as the "tag side"), then look for consistency with an 
invisible decay or a decay to a single photon of the other neutral B (referred to as 
the "signal side"). 

In the signal event selection we consider events with no charged tracks besides 
those from the B — > D^*> + £~V candidate. In order to reject background events 
where one charged or neutral particle is lost along the beam pipe, the cosine of the 
polar angle of the missing momentum in the CM frame (cos0^ ss ) is required to lie in 
the [—0.9, 0.9] range. For the B° — > invisible decay, in events where the D meson on 
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Figure 2: Distributions of the NN output for simulated B° — > invisible events with a 
D meson on the tag side. The black solid line is the signal while the red dashed line 
is the background. The solid gray vertical line defines the NN output signal region. 

the tag side decays into K~7r + -K~, two additional selection criteria are also applied: 
on the sum of the angles between the Kaon and each one of the two its and on the 
sum of the angles between the lepton and each one of the two 7rs. To reconstruct 
B° invisible+7 events, one remaining photon candidate with energy greater than 
1.2 GeV in the CM frame is also required. 
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An artificial neural network (NN) is used to provide further discrimination between 
signal and background events. Events with a D or a D* meson on the tag side are 
split in two different categories and a different NN was used for each sub-sample. 
For the B° — y invisible decay 9 and 6 variables are used for D and D* sub-samples, 
respectively, while for the B° —y invisible+7 decay we used 6 and 4 variables. These 
variables are mainly kinematical and refer to the tag side reconstruction. The two 
most important are the cosine of the angle between the B meson and the D { *H pair, 
defined as 

2E B E DM - e+ - m| - m 2 D(t) _ g+ 

cose BDM - t+ = — y — s (1) 

2\PB\\PD(.*)-i+\ 

and M^iss (defined as the invariant mass of the event after the D^~£ + pair is sub- 
tracted). In the B° — y invisible+7 analysis, we additionally use the energy of the 
photon on the signal side, evaluated in the laboratory frame. In Fig. [2J the output 
of the NN for simulated B° — y invisible with a D meson on the tag side, and the 
corresponding signal region, are shown. 

After the NN selection, the D meson invariant mass (itid) and the difference 
between the reconstructed D* invariant mass and the PDG D° mass (Am) are used 
to define a signal region and a sideband region for the D tag and D* tag samples, 




Figure 3: Results of the maximum likelihood fit of £" e xtra for B° —y invisible (left) and 
B° —y invisible+7 (right). 

respectively. The signal region is defined as a ±15 MeV/c 2 window around the PDG 
value for mo, and as 0.139 < Am < 0.148 GeV/c 2 . The excluded regions are used as 
sidebands. 
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The neutral energy that remains after all tag side tracks and neutral clusters have 
been accounted for is denoted as -E ex tra- For B° — > invisible+7, the energy of the 
highest-energy photon remaining in the event (the signal photon candidate) is also 
removed from the -E ex tra computation. The -E ex t m signal region is defined by imposing 
an upper bound at 1.2 GeV. 

We construct probability density functions (PDFs) for the E extia distribution for 
signal (P S ig) and background (Pbkg) using detailed MC simulations for signal and 
data from the mo and Am sidebands for background. The two PDFs are combined 
into an extended maximum likelihood function C, defined as a function of the free 
parameters iV s i g and A^kg, the number of signal and background events, respectively. 
The photon reconstruction algorithm has a lower cluster energy cut of 30 MeV, and 
as a consequence, the -E ex tra distribution is not continuous. To account for this effect, 
the likelihood is composed of two distinct parts, one for E extT & > 30 MeV and one for 
£ ex tra = MeV. 

The negative log-likelihood is then minimized with respect to N sig and A^kg in 
the data sample. The resulting fitted values for N sig and A^kg are given in Table [TJ 
Figure [3] shows the -E ex tra distributions for B° — > invisible and B° — > invisible+7 with 
the fit superimposed. 



Mode 


N ■ 


AW 


B° — > invisible 


-22 + 9 


334 ± 21 


B° — > invisible+7 


-3.1+5.2 


113 + 12 



Table 1: Fitted yields of signal and background events in data. The uncertainties 
are statistical. 

Using detailed Monte Carlo simulations of B° — > invisible and B° — > invisible+7 
events, we determine our signal efficiency to be (17.8 ± 0.2) x 10~ 4 for B° — > invisible 
and (16.0 ± 0.2) x 10~ 4 for B° — > invisible+7, where the uncertainties are statistical. 



Source 


B° ->• invisible B" - 


-> invisible+7 


Source 


B" — > invisible 


B° — > invisible+7 


Normalization Errors 




Efficiency Errors 




B-counting 


0.6% 


0.6% 


Tagging Efficiency 


3.5% 


3.5% 


Yield Errors (events) 




mn (Am) Selection 


1% 


1.3% 


Background Param. 


15.8 


6.5 


Preselection 


3% 


2.4% 


Signal Param. 


2.0 


1.2 


Neural Network 


6.1% 


8.2% 


Fit Technique 




1.0 


Single Photon 




1.8% 


-Eextra Shape 


0.1 


1.8 


TOTAL 


7.7% 


9.5% 


TOTAL 


15.9 


6.9 





Table 2: Summary of the systematic uncertainties. 



The systematic uncertainty on the signal efficiency is dominated by data-MC dis- 
crepancies in the distribution of the variables used as input to the NN while the 
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main systematic uncertainty on the signal yield is dominated by the background 
parametrization uncertainties. The total systematic uncertainty on the signal selec- 
tion efficiency is 7.7% for B° — > invisible decay and 9.5% for B° — » invisible+7 decay 
and the total systematic errors on the signal yield are 16 and 7 events for B° — > invis- 
ible and B° — > invisible+7, respectively. All systematic uncertainties are summarized 
in Table H 

3 Conclusions 

A Bayesian approach is used to set 90% confidence level (CL) upper limits on the 
branching fractions for B° — > invisible and B° — > invisible+7. Flat prior probabilities 
are assumed for positive values of both branching fractions. Gaussian likelihoods are 
adopted for signal yields. The Gaussian widths are fixed to the sum in quadrature of 
the statistical and systematic yield errors. We extract a posterior PDF using Bayes' 
theorem, including in the calculation the effect of systematic uncertainties associated 
with the efficiencies and the normalizations, modeled by Gaussian PDFs. Given the 
observed yields in Table [[J the 90% CL upper limits are 

B(B° invisible) < 2.4 x 10" 5 
B(B° invisible + 7) < 1.7 x 10" 5 

at 90% CL. These limits supercede our earlier results [8], which used a small fraction 
of our present dataset, and the recent Belle result [9]. 

References 

[1] S. Adler et al. (BNL-787 Collaboration), Phys. Rev. Lett. 84, 3768 (2007); 

[2] A. Dedes, H. Dreiner, and P. Richardson, Phys. Rev. D 65, 015001 (2002). 

[3] K. Agashe, N. G. Deshpande, and G.-H. Wu, Phys. Lett. B 489, 367 (2000). 

[4] K. Agashe and G.-H. Wu, Phys. Lett. B 498, 230 (2001). 

[5] H. Davoudiasl, P. Langacker, and M. Perelstein, Phys. Rev. D 65, 105015 (2002). 

[6] B. Aubert et al. (BABAR Collaboration), Nucl. Instrum. Meth. A 479 (2002). 

[7] J. P. Lees et al. (BABAR Collaboration), Phys. Rev. D 86, 051105 (2012). 

[8] B. Aubert et al. (BABARCollaboration), Phys. Rev. Lett. 93, 091802 (2004); 

[9] C. L. Hsu et al. [Belle Collaboration], Phys. Rev. D 86, 032002 (2012). 



5 



